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Abstract. The decay of photoexcited quantum systems (examples are 
photodissociation of molecules and autoionization of atoms) can be viewed as a 
half-collision process (an incoming photon excites the system which subsequently 
decays by dissociation or autoionization). For this reason, the standard 
statistical approach to quantum scattering, originally developed to describe 
nuclear compound reactions, is not directly applicable. Using an alternative 
approach, correlations and fluctuations of obscrvablcs characterizing this process 
were first derived in [Fyodorov YV and Alhassid Y 1998 Phys. Rev. A 58 R3375]. 
Here we show how the results cited above, and more recent results incorporating 
direct decay processes, can be obtained from the standard statistical scattering 
approach by introducing one additional channel. 
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Molecular photodissociation ^ [21 01 01 [HI EH and atomic autoionization [7] are 
examples of quantum-mechanical decay processes: excited quantum systems decay 
into one or several scattering continua (decay channels). The situation is illustrated 
schematically in figure ^ An atom or a molecule is excited from its ground state or 
a low lying state \g) by a photon. Photoabsorption cither results in direct excitation 
of the decay channels \E a ) (where a = 1,...,M) or in excitation of quasi-bound 
states |n) (where n = 1, . . . , N). Decay into a channel \E a ) then proceeds via residual 
interactions between the quasi-bound states and the continua. 

Correlations and fluctuations of the quasi-bound states determine the statistical 
properties of the indirect decay cross section. If the corresponding classical 
dynamics is sufficiently chaotic, these properties can be modeled by random matrix 
theory [HI El I1C)| . In reference this approach was adopted to characterize the 
distribution and two-point correlations of the indirect photodissociation cross-section, 
in the absence of direct decay channels. The results of m were obtained by calculating 
the statistical properties of a resolvent of an effective non-Hcrmitian Hamiltonian. 
Equivalently they may be obtained from eigenvector correlations |12l 1131 IT^] of the 
non-Hcrmitian effective Hamiltonian. In the limit of weak continuum coupling, cross- 
section correlations can be calculated using perturbation theory [Uj. The results of 
reference |15] have been generalized in |16j where interference effects between direct 
and indirect decay paths in the photoabsorption cross section were studied. 

How are the results obtained in |1 HI15IITH] related to those derived for nuclear and 
mesoscopic scattering processes within the so-called statistical scattering model (also 
termed the Heidelberg model ^1 EH1 E31)? Since photodecay is a half-collision 
process [201 l2~T] . the statistical scattering model as such is not directly applicable. 
In order to establish a connection to the statistical scattering model, it has been 
suggested [221 to introduce an additional channel representing photoexcitation. In 
this paper we show how to describe the decay of photoexcited quantum systems using 
this idea. Our results can be summarized as follows: We demonstrate that partial 
and total decay cross sections can be obtained from an extended (almost) unitary 
scattering matrix by introducing an additional channel for photoexcitation and a linear 
transformation reminiscent of the "Engelbrecht-Weidenmuller" transformation 
We use that relation to obtain a non-perturbative result for the two-point correlation 
function of the total cross section in a random-matrix model, which allows for arbitrary 
coupling strengths between quasi-bound states, decay channels and the photoexcited 
initial state. 

Let | a) = jl\g) be the ground state dressed with the dipole operator jl. The 
partial photodissociation (photoionization) cross section into channel a is proportional 



where the scattering wave function |\Ja°)) is outgoing in channel a only. The total 
photodissociation cross section is given by 

M 



(M is the number of open channels). We follow the approach adopted in ^2 DS| 
where resonant and direct processes were treated separately by means of the Feshbach 
approach |2S]: the Hilbert space of excited states is divided into a subspace of 



to 122121 



(E)<x (a\^ a \E)) 



(1) 




(2) 



a=l 
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Figure 1. Direct photoionization from an initial state \g) into the subspace 
of unbound states (M continuum channels \E a ), here: M = 3) with transition 
amplitudes a° ut , and indirect photoionization via N bound states \n) with 
transition amplitudes aJJ 1 . The states \n) turn into resonances due to couplings 
to the continua with matrix elements W na - 



bound states, and its complement containing all unbound states. The vector \a) is 
decomposed accordingly | a) = |a m ) © |a out ). Choosing the bound space appropriately, 
the scattering problem projected onto the complement looses all resonant features. 
The Hamilton operator in the bound space is represented by an N x N random matrix 
Ho, with eigenstates \n), n = 1, . . . N. The N dipole transition matrix elements to 



are denoted by a™, and a™ = (a 1 " 



<N 



) J . Correspondingly, the M (energy- 



independent) dipole transition matrix elements from \g) to the decay channels are 
a° ut , and cc ' 



„out\T 



i M ) . The matrix elements of the residual interaction 
W are denoted by W na = (n\W\E a ), and W a = (Wia, ■ ■ ■ , Wno,) t arc the column 
vectors of W . For the total cross section, the following result has been obtained in 



OS: 



a{E) = ||a out f -TT-Hm 



(a in + i7r Wa out )i 



1 



(a in -mW a out ) 



(3) 



E — H c ff 

Here H c h = Hq — inWW'' is an effective non-Hcrmitian Hamiltonian. The partial 
cross section is given by 

2 



aa(E) 



wl 



1 



a E H, 



oil' 



(a i: 



in W a c 



(4) 



It can be shown by purely algebraic manipulations that, indeed, X^^li a a{E) = <r(e). 

We now prove that the cross sections given in J3Jl and can be expressed in 
terms of the extended (1 + M) x (1 + M) unitary scattering matrix 



S(E,S) = l 1+M -2m V f 



E-R 



efi 



V with V = (cc in 5/2tt , W) , 



(5) 



Letter to the Editor 



4 



and F e g = Hq — inV V' . This maps our problem onto the statistical scattering 
model. It is convenient to enumerate the channels from to M . so that the new 
channel, referring to photoexcitation, gets the index 0. The corresponding column 
vector of V is scaled with a real positive parameter 8 (which will be taken to zero at 
the end of the calculation). Now, consider transform S'(E,S) = u(5) S (E , 8) v T (8) , 
where the matrices u(8) and v(8) are given by 



u(8) 



i 


(a out )t 8/2 \ 




( 1 


(a out ) T 8/2 \ 







v(8) = 





1a/ 





J 




^ o 





(6) 



The transformation © is reminiscent of an "Engclbrccht-Wcidcnmullcr" transforma- 
tion although for <5 > 0, u(8) and v(6) are not unitary. We now show that the 
cross sections © and @ can be expressed in terms of the elements of S'(E,S) as 
follows 

2 



0- a (E) 

a(E) 



lim <r 2 

2 lim <T 2 

5^0 



-a^ + S' a0 {E,5) 



1 + -RbS' 00 (E, 5) 



(7) 

(8) 



In order to derive these relations, we compute the corresponding matrix elements of 
S'(E,8). Let 1 < a < M: 

M M 



SUE, S) 



u ab (8) S bc (E, 8) v 0c (S) = i S a0 (E, 8)+J2 S ac (E, 8) - a° 

b,c=0 c=l 
S out + 1 



8 W, 



a E-F cS 

c A/ 



(a m -inW a out ) 



(9) 



S' m (E, 8) = -Soo(E, S)+i-J2 [(<0* S ba (E, 8) + S ab (E, 8) <*] 



X2 M 
6,c=l 

8 2 



b=l 

l^)*S bc {E,8) a° 
1 



-1 + i _ (a in )t 

2tt 1 ' E-F eS 



a 



1 



a m + (a m ) 



1 



~ 1 + T " a ° 



E — F e s E — F e g 

i 8 2 1 



(10) 



2tt ' E 



oir 



(a> - Itt W a out ) 



Since i^ff — > i? c ff in the limit (5 — > 0, we obtain expressions © and 10} for the cross 
sections. This is easily verified by inserting JjJJ and (|10fl into the left hand side of Q 
and respectively. 

Equations J7J and (JSJ) enable us to compute the average total cross section, as 
well as its correlation function, defined as 

C(E, w) = (a(E - w A/2) a{E + w A/2)} - {a{E)) 2 . (11) 

Here A is the average local distance between neighboring resonances, and (. . .) denotes 
an energy average, often replaced by an ensemble average in theory For simplicity, 
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we assume that the column vectors of V in equation (jSJ) are pairwisc orthogonal (see 
below). In this case, the average extended scattering matrix (S(E,5)) is real and 
diagonal, with diagonal elements given by [5^1 

{S aa (E,5)) = l-^, A o =7r 2 0o||W (a) || 2 , and A = ^ go ||a in || 2 , (12) 



1 + A„ 

where Nqq = 1/A is the level density. Starting from equation (JHJ) and assuming 
that Hq corresponds to a classically chaotic system, we obtain the expression for the 
average total cross section derived in ^Hj : 

M 

1 + ^ Re U Qa V 0b (S a b(E, 8)) 



HE)) 



\a 



lim — 

s^o S 2 



a,b=0 



A I 



= a 



out||2 



Po||« in || 2 -E 



„out 1 2 



A„ 



1 + A a 



(13) 



In a similar way, we may write for the autocorrelation function of o~(E) 

2 M M 

C(E,w) = \im — ^2 ^2 u 0a(S) u 0c (S) v ob (5) v 0d (S) 



x Re 



[(S ab (E - wA/2) S* cd (E + wA/2)) - (S ab (E)) (S* cd (E))] . (14) 



It has been used that the correlation function of two 5-matrix elements is different 
from zero only if one S'-matrix element is complex conjugated while the other is not 
(see reference j^Hl for details). Equation ifH|) expresses the autocorrelation function 
of the total absorption cross section in terms of correlation functions of elements of 
the extended scattering matrix, defined in J^J). Usually, the energy dependence of 
the autocorrelation function is mainly due to the energy dependence of the coupling 
parameters a out , ||o! m || 2 and Ai, . . . , Am- This dependence should be weak in order 
to allow energy averaging over a sufficiently large window. In cases where the 
corresponding S-matrix correlation functions are available, equation (|14fl enables us 
to obtain the autocorrelation function of the total cross section in closed form. 

As an illustration, let us assume time-reversal invariance for the collisional system 
and predominantly chaotic dynamics in the large (N —> oo) bound space. Then, the 
required S'-matrix correlation functions are given by the Verbaarschot-Wcidcnmullcr- 
Zirnbaucr (VWZ) integral [18) ■ Using a Fourier representation developed in we 
obtain: 



C(E,w)=4 J At e- 27ri wt J dr J 



max(0,i — 1) 



du , 
o (2« 



(t - r)(r + 1 - t) 



l)(t* 



M 

n 



1 - T e (t - r) 
\ (1 + 2^ + ^)1/2 




M 



n 



nil 



M 



TgT b 

16 



n„ 



(15) 



where x = u 2 (2r + l)/(2u + 1), and where r a = T a \a a ut \ 2 /(g ||a in || 2 ). As usual, 
it is more convenient to work with the transmission coefficients T a = 4A Q /(1 + A a ) 2 
to describe the coupling between the subspace of quasi-bound states and the decay 
channels. The functions A a and Tl a b arc defined as follows: 
r + T„ x t — r 



1 + T a (2r + T a x) l-T a (t-r) 
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n 



ab 



T a T b x 1 + [T a T b r+(T a + T b )(r + 1) - l]x + r(2r + 1) 



(1 



2T Q r 



r)(r + 1 



2T b r 
-t) 



T b x) 



(16) 



■ r)][l-T b {t-r)] ■ 
1 open channel, a result corresponding to (15,16) was 



[1 - T a (t 

In the special case of M - 
derived, by a different approach, in |27j . In the absence of direct processes (i.e., when 
all direct dipolc transition amplitudes a° ut are zero) equation (15) reduces to the 
result derived in JT] . Note that in l|16|) , the limit S — > has been taken. This implies 
that the transmission coefficient corresponding to the "photo channel" has to be set 
to zero. 

In figure the correlation function C(E,u>) is shown: as obtained from 
equation (15), valid in the limit of N — > oo, and as obtained from the numerical 
diagonalisation of finite matrices of size N = 128. In the presence of direct decay 
channels, finite-size effects are seen to be somewhat larger, as compared with purely 
indirect decay. 

We conclude with two remarks. First, the VWZ-integral can also be used to 
calculate the average of the partial photoabsorption cross sections. However, for 
the correlation functions of the partial cross sections no exact analytical formulae 
are known. In this case one must use approximate results, such as the "rescalcd" 
Breit-Wigner approximation |1 51 126| . Second, in theoretical studies of the statistical 
scattering model, it is often assumed that the channel vectors (the column vectors 
of V) arc orthogonal to each other: one can always obtain the original scattering 
matrix by an appropriate unitary transformation from a simpler scattering matrix 
with orthogonal channel vectors. 
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Figure 2. Correlation function 1111 of the total cross section, time-reversal 
invariant case, results according to equation (15) (lines) and from the numerical 
diagonalisation of random matrices with N = 128 (symbols), (a): without direct 
coupling (c* out = 0), M = 2, 7 = 0.02 (o), and 7 = 0.5 (□). (b): with direct 
coupling, M = 2, 7 = 0.02, c^ ut = 2.5 (o), and 7 = 0.5, a° ut = 0.5 (□). 
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